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We study the penguin pollution in the B° — > J/ipKs decay up to leading power in 1/rrib and 
to next-to- leading order in a s , mi, being the 6 quark mass and a 3 the strong coupling constant. 
The deviation ASj/^,k s of the mixing-induced CP asymmetry from sin(20i) and the direct CP 
asymmetry Aj/^k s are both found to be of O(10 -3 ) in a formalism that combines the QCD- 
improved factorization and perturbative QCD approaches. 



I. INTRODUCTION 

The B° — > J/tpKs decay is known to be the golden 
mode for extracting sin(20!), 4>i being the weak phase of 
the Cabbibo-Kobayashi-Maskawa (CKM) matrix element 
Vtd, through the time-dependent CP asymmetry 



a(t) = 



T(B°(t) -> J/^K S ) - T(B°(t) -> J/^K S ) 
r(S°(t) -> J/ipKs) + T{B°(t) J/^K S ) ' 
S J/ilKs sin(AMi) + A,,^ Ks cos(AAft) , (1) 



where AM is the mass difference of the two i^-meson 
mass eigenstates. The mixing-induced CP asymmetry 
Sj/^Ks i s naively equal to sin(20i) in the standard 
model. At the same time, the direct CP asymmetry 
Aj/^Ks is expected to be vanishingly small. The pen- 
guin pollution, which would change the above naive ex- 
pectations, is believed to be negligible in this mode. A 
complete estimation of its effect is, however, essential for 
future precision measurements, where the experimental 
error of Sji^k s will be 5Sj/^x s ~ ±0.2 for 5ab _1 data 
at the Super B factory 0. 



The deviation AS 



'J/TpKs 



sin(20i 



found to be of 0(1CP 4 ) in the previous estimation by 
Boos et al. Q, taking into account corrections to the 
B — B mixing and to the B° — ► J/ipK° decay amplitude. 
A part of penguin-operator contributions, however, were 
overlooked in the estimation. 

A model-independent approach to estimate the size of 
the penguin pollution proposed by Ciuchini et al. 
They extracted the range of the penguin corrections from 
the data of the branching ratio and the CP asymmetries 
for the B° — > J/ipTr° decay, and used the range to evalu- 



ate AS 



J/^Ks 



with the flavor SU(3) symmetry. The cur- 



rent data lead to 



AS j C /7k s = 0-000 ±0.012 , 



(2) 



where the large error comes from the experimental un- 
certainty of the B° — ► J I ipn data. 




FIG. 1: Non-local corrections to the B- 
dots represent local operators. 



B mixing, where the 



In this talk, we present the most complete analysis 
of the penguin pollution in B° —> J '/ipKs up to leading 
power in l/m& and to next-to- leading order (NLO) in a s , 
nib being the b quark mass and a s the strong coupling 
constant. We shall concentrate only on the corrections to 
the B — B mixing and to the decay amplitude, and refer 
the inclusion of corrections to the K — K mixing and to 
the _B-decay- width difference Ar^ to Ref. 

This talk is organized as follows: In Sec. HH we briefly 
review the corrections to the B — B mixing. The cor- 
rections to the decay amplitude are studied in Sec. IIIII 
Numerical results are then presented in Sec. IIVI Sec- 
tion W\ is the conclusion. 



II. CORRECTIONS TO THE B — B MIXING 

The mixing-induced CP asymmetry Sj/^k s an d the 
direct CP asymmetry Aj/^ Ks in Eq. |T]) are given by 



S 



21mA 



J/ipKs 



J/4>K S 



1 + \^J/il>Ks\ 
\Kl/jpKs\_ ~ 1 
1 + \^J/^Ks\ 2 



with the associated factor 

q A(B° -> J/ipK S ) 



A 



J/ipK s 



p A{B° -> J/4>K S 



(3) 



(4) 
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where the ratio q/p is related to the B — B mixing and 
A(B°(B°) -> J/ipK_ s ) is the decay amplitude. The cor- 
rections to the B — B mixing alters q/p from exp(— 2i(f>i). 
In Ref. [1} , non-local contributions to the B — B mixing, 
shown in Fig. [TJ were calculated. The non-local operators 



are defined as 



^xT[[{bc)^cdY]{x)[{lc) v {cdY] 



T-2 = --J d 4 xT [[(bc)^cdr](x)[(bu) v (udy](0)] , 
T 3 = -- [ d 4 xT [l(bu)^udr](x)[(buUudy](0)] , (5) 

where the left-handed current is shortened into {bq)^ = 
(Pl1h<Il), and T denotes the time-ordered product. As 
the scale evolves down, the non-local operators mix into 
the AB — 2 local operators of dimension 8, 



Qr = D(bd)^bdr , 

q 2 = d»d lJ {bd)^bd) 1J 

Q 3 = m 2 c {bd) tl {bdy , 



(6) 



m c being the c-quark mass, through renormalization- 
group effects. Including QCD corrections at 0(a s ) and 
resumming the large logarithm a s hi(mw /irib), where 
mw is the IF-boson mass, the corrections to the mix- 
ing factor q/p lead to [|[ 



a nmix 

a mix 

A J/ipK s 



(2.08 ± 1.23) x 1(T 4 
(2.59 ± 1.48) x 1(T 4 



(7) 



Thus, AS™}* Ks and A v ^ Kg have been both found to be 
ofO(l(T 4 ). 



III. CORRECTIONS TO THE DECAY 
AMPLITUDE 



The B° — > J/ipK° decay amplitude can be decom- 
posed into 

A(B° -> J/ipK°) 

T/* T/ /|( u ) _i_T/*t/ /|( c ) i T/*T/ 

— v ub v us ^j/^ko "r r cfc •'es -^j/^ko t" I'tfa I'ts -^j/^ifO , 



(*) 



(8) 



where the unitarity relation of the CKM matrix ele- 
ments V t * b V ts = -V* b V us - V* b V cs is used. The term 
with V* b V cs dominates the decay amplitude, whereas 
that with V* b V us causes the penguin pollution, since it 
carries the CKM phase 3 , which is defined by V ub = 
\V u b\ exp(-i(f> 3 ). Note that the amplitude V* b V us A ( jj^ K0 
was missed in Ref. 

We calculate the decay amplitude using a formal- 
ism Q that combines the QCD-improved factoriza- 
tion (QCDF) 7] and the perturbative QCD approach 
(PQCD) [8( up to leading power in l/m b and to NLO 
in a s . We adopt QCDF to handle the factorizablc 
amplitudes, since the energy release in the form factor 
F+ K (m 2 ^), mj/^p being the J/tp meson mass, is small 
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FIG. 2: Factorizable diagrams, where the spectator quark is 
omitted. 



and it is unlikely to further factorize the form factor. 



The factorizable contributions to A 



(c) 

J/ipK° 



and A 



(*) 

J/ijjK 



displayed in Fig. [5] are given by 
Af^ K0 = 2V&J dx c * L (x c )F* K {m 2 ^)a 2 {x c ,t) , 



■^j/Iko = 2V6 [ dx c ^ L (x c )F^ K (m 2 Jfi ,) 
x [a 3 {x c ,t) + a 5 (x c ,t)] 



(9) 



respectively, where x c is the momentum fraction carried 
by the c quark, and 1 i' L (x c ) the twist-2 J/tf) meson dis- 
tribution amplitude. The explicit expressions of i Sf L (x c ) 
is referred to Ref. [l|. The hard scale t has been chosen 
as 



t = max 



\\J x o{m B -m^), ^x c (m 2 B ~m 2 ^) ) , (10) 

with the notation x^ = 1 — x c and tub being the B me- 
son mass. The effective Wilson coefficients, including the 
0(a s ) vertex corrections in Figs. [5£b)-(d), are defined as 



a 2 (x,fi) = Cx(fx) + 



Nr. 



+ ^M CF (_ 18 + 12]n ^ + fl(x) 

a s (x,(i) = CaOi) + C 9 {iM) + — [Ci{/i) + Ci (jM)] 

+ ^C F (-18 + 12^+f I (x) 
a 5 (x,n) = C 5 (fi) + C r (fJ,) + — [Ce(fx) + C s (fi)] 



1 H -. — C F 6 - 12 In fi(x) 

An \ [i 



(11) 



with N c ~ 3 being the number of colors, Cp = 4/3 the 
color factor, and Ci_io the Wilson coefficients defined, 
e.g., in Ref. @. The loop function fi(x) is given by @ 



//(*) 



3(1 - 2x) 



1-x 



In x — 37ri 



+31n(l - r 2 



1 - r J/*P X 



■(12) 





d d d d 

FIG. 3: Nonfactorizable spectator diagrams. 





FIG. 4: it-quark loop diagrams (a) with a photon emission 
and (b) with a gluon emission. 



with the ratio rj/^ = mj/^/rriB- 

For the 0(a s ) nonfactorizable spectator diagrams in 
Fig. [3l QCDF is not appropriate due to the end-point 
singularity from vanishing parton momenta. Note that 
the nonfactorizable contribution has a characteristic scale 
higher than in {rn?j,^) [101 ] . Therefore, we can em- 
ploy PQCD based on kx factorization theorem, which 
is free of the end-point singularity. The nonfactorizable 
spectator amplitudes in PQCD are given by 



/)(<=)"/ 

At)nf 



= M (4) , 
= M(a' 3 -a' 5 ) 



(13) 



where the explicit expression of the amplitude M. is re- 
ferred to Refs. p], @ , and the effective Wilson coefficients 
are defined by 



40) 
40) 
40) 



Cap) 



1 

X 
1 

N~c 



[c 4 0) + CioO)] 
[c 6 0) + c 8 0)]- 



(14) 



The nonfactorizable contribution is essential for the 
B° — > J/%jjK° decay, since this mode is classified in the 
color-suppressed category of B meson decays. In fact, 
the nonfactorizable contribution is comparable in size to 
the factorizable one. 

The amplitude -^j^ K o receives the contribution from 
the u-quarks loop shown in Figs. Ufa) and (b). The 
penguin corrections from the it-quark loops as well as 
from the c-quark ones are well-behaved in perturbation 
theory without any infrared singularity, which has been 
known as the Bander-Silverman-Soni mechanism (Tlj . It 
means that their long-distance contribution is unlikely 
to be large. The light-cone-sum-rules analysis has also 




(a) 



(b) 



FIG. 5: QCD tt-quark loop diagrams (a) with an additional 
hard gluon and (b) with a soft gluon. 



suggested that these corrections are dominated by short- 
distance contribution The u-quark loop contribu- 
tion was estimated naively as ASj° c ^ (s = —(4.24 ± 

1.94) x 10" 4 0. We shall reinvestigate this contribution 
below. 



The it-quark loop contribution can be expressed in 
terms of the effective Hamiltonian 



n 



(«) 

efi 



3tt 



-j=V us V ub I - In • 



e u e c (iVcdO) + C 2 (/x)) (c 7ai c) (s 7 "(l - 75 )6) 



o., 



+ ^O) (c lfi T a c) - l5 )T%) 



(15) 



where Gp is the Fermi constant, a the fine-structure con- 
stant, e n ( c ) — 2/3 the electric charge of the u(c) quark, 
and T a the SU(3) generator. The first and second terms 
arise from the photon emission diagram in Fig. IHa) an< ^ 
the gluon emission one in Fig. |4] (b), respectively. I 2 
denotes the invariant mass squared of the photon or 
gluon. The photon emission contribution is less than 5% 
of A < 'jJ^ K due to the smallness of a, and can be safely 
dropped. For the gluon emission, an additional gluon is 
necessary to form the color-singlet J/ip meson. If the 
additional gluon is hard, the resultant contribution in 
Fig. [UJa) is of next-to-next-to-leading order in a s and 
beyond the scope of this work. If the additional gluon 
is soft, the corresponding nonperturbative input is the 
three-parton J/ ijj meson distribution amplitude as shown 
in Fig. [5^b). The twist-3 distribution amplitude, which 
contributes only to transversely polarized J/tj) mesons, 
is irrelevant here. The twist-4 distribution amplitudes 
are antisymmetric under the exchange of the momentum 
fractions of the c quark and of the c quark. Because the 
hard kernel associated with the it-quark loop is symmet- 
ric under the above exchange, its convolution with the 
distribution amplitudes diminishes. Therefore, the gluon 
emission contribution as well as the photon emission one 
can be neglected. The c-quark loop corrections can also 
be neglected, since it modifies only the branching ratio 
slightly. 



4 



In summary, the amplitudes in Eq. {8j are given by 
, 



A^ 



0) 

J/ij>K a 
J/ipK° 



J/i/jK° 

(*)/ 
J/ipK° 



A 
A 



{c)nf 
J/ipK° 

(t)nf 
J/i>K a 



(16) 



where A {c ^ f and A {c ' t]nf are found in Eqs. © and CL1, 
respectively. 



Including the correction to the B—B mixing in Eq. ([7]), 
we predict ASj/^Ks an d ^-J/tpK s ^ U P to leading-power in 
and to NLO in a s , as 



AS 1 



i J/tpK s 



a cmix _,_ a odccay 



(9.3+ 36 



4.6 



-4 



;) xio 

* decay 
" ^J/^Ks 



,4 mix 

A J/ipK s 

-(I4.1^ 2 )xl0- 4 . 



(18) 



IV. NUMERICAL RESULTS 

We predict the branching ratio, ASj/^k S i an d 
Aj/ipK S i with the corrections to the decay amplitude, 



Br( J B° -> J/^ ) 



6.6 



AS" 



decay 
J/ipKs 



7.2 



+3.7 (+3.7) 
-2.3 (-2.3) 

+2.4 (+1.2) 
-3.4(-l.l) 



x 10- 4 , 
x 10- 4 , 



* decay _ 
A Jli>K s ~ 



16.7 



+6.6 (+3.8) 
-8.7(-4.1) 



X 10" 



(17) 



where the choices of parameters used in the computation, 
including their allowed ranges, are referred to Ref. [l[. 
The errors are estimated from the variation of the CKM 
and hadronic parameters, while those in the parenthe- 
ses arise only from the hadronic ones. The predicted 
branching ratio is in agreement with the data Br(S° — > 
J/^K°) = (8.63 ± 0.35) x 10" 4 [13]. The form factor 
F+ K (m 2 jj^) increased by 15% and larger Gegenbauer 
moments in the kaon distribution amplitudes can easily 
account for the central value of the data. Our ASj°^ s 

( A jJ^k s ) from the penguin corrections is about twice of 
the naive estimation from the w-quark loops [|[ , and has 
an opposite (the same) sign. Both AS"jy^ s and A^° c ^ (s 

in Eq. (fT7|) indicate the O(10~ 3 ) penguin pollution in the 
B° — > JhhKs decay, consistent with the conjecture made 
in Ref. 



Taking into account the CP violation from the K — K 
mixing Q, ASj/^k s and Aj/^, Ks remain O(10~ 3 ). 



V. CONCLUSION 

In this talk, we have presented the most complete anal- 
ysis of the branching ratio and the CP asymmetries for 
the B° — > J/ipKs decay up to leading power in 1 jm^ and 



to NLO in a s . We have found that both AS 



J/ipKs 



and 



.4 



are of O(10 ). Our prediction for AS, 



is smaller than an expected systematic error in future 
data As for the direct CP asymmetry, our result 

supports the claim that Aj/^ Ks > 1 % would indicate 
new physics [lj| • These predictions provide an important 
standard-model reference for verifying new physics from 
the B° -> J/tpKs data. 
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